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Abstract. This paper deals with flow-induced shape transitions of elastic capsules. 
The state of the art concerning both theory and experiments is briefly reviewed starting 
with dynamically induced small deformation of initially spherical capsules and the 
formation of wrinkles on polymerized membranes. Initially non-spherical capsules 
show tumbling and tank-treading motion in shear flow. Theoretical descriptions of the 
transition between these two types of motion assuming a fixed shape are at variance 
with the full capsule dynamics obtained numerically. To resolve the discrepancy, we 
expand the exact equations of motion for small deformations and find that shape 
changes play a dominant role. We classify the dynamical phase transitions and obtain 
numerical and analytical results for the phase boundaries as a function of viscosity 
contrast, shear and elongational flow rate. We conclude with perspectives on time- 
dependent flow, on shear-induced unbinding from surfaces, on the role of thermal 
fluctuations, and on applying the concepts of stochastic thermodynamics to these 
systems. 

PACS numbers: 47.15.G- (Low-Reynolds-number (creeping) flows), 87.16. D- 
(Membranes, bilayers, and vesicles), 83.50.-v (Deformation and flow) 
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1. Introduction 

Elastic capsules constitute an important class within soft matter research both from 
a fundamental and application point of view |TJ [2j [3j HJ E]. In contrast to strongly 
deformable fluid vesicles [B], the membrane of capsules exhibits finite shear elasticity. 
This structural feature is achieved either by chemically or physically crosslinking the 
molecules forming the membrane [7j or by attaching a polymeric network to a fluid 
membrane as in red blood cells (RBCs) [S]. From the perspective of this paper 
these biological capsules can be subsumed under this heading of a micro-capsule as 
well. Promising applications of these objects arise interalia from the pharmaceutical 
perspective of optimized drug encapsulation delivery and release [91 IIP]. 

In a microfluidic context, one of the most relevant issues is the behaviour of 
these soft objects in externally applied flow. This problem is challenging from a 
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theoretical perspective since the shape of the capsule is not given a priori but determined 
dynamically from a balance of interfacial forces with fluid stresses. Hence, their 
deformation will depend on their material properties and the type of external flow. 

Before reviewing the state of the art concerning theory and experiment in the 
next two sections, we set the stage by introducing a few key notions for micro-capsule 
dynamics. As illustrated in figure [TJ we consider the dynamics of a single closed elastic 
membrane that surrounds an inner fluid with viscosity r)i and is embedded in an ambient 
flow with viscosity r] , thereby defining an important quantity, the viscosity contrast 
i]i/r) . In the absence of the capsule we assume a prescribed external flow, in most cases 
set to be an infinite linear shear flow. For all experimental setups the Reynolds number 
is small, i.e. dissipation is dominating and inertial effects are negligible. Consequently, 
the flow is governed by the linear Stokes equations and determined instantaneously by 
the boundary conditions, where in all cases a no-slip boundary condition is assumed at 
the membrane. 

For describing the capsule state adequately, we define a set of suitable parameters. 
Capsules and RBCs consist of a thin membrane with a typical thickness in the nanometer 
range which encloses a fluid on a larger length scale (typically of the order 1-100 [im) 
with fixed volume V defining a length scale R as the radius of a sphere with the same 
volume 

V = grf. (1) 

For capsules with area elasticity, the area A of the membrane, as well as the excess area 
A > defined by 

A = (47r + A)R 2 , (2) 

can change, where A = corresponds to a sphere. In contrast, the membrane of RBCs 
and polymerized vesicles has a vanishingly small area compressibility, i.e. A and A are 
constants. For these area-incompressible objects, shape changes are only possible for 
initially non-spherical shapes A > 0. 

In most cases considered in this paper, the capsule is roughly ellipsoidal, see figure [T] 
for a definition of the geometry. The three main dynamical parameters are (i) the 
inclination angle \1/ between the direction of the shear flow and the long semi-axis, (ii) 
the phase angle $ measuring the angle between the direction of the shear flow and 
that of the position of a tracer particle which initially was on the endpoint of the long 
semi-axis and (iii) the Taylor deformation parameter 

D = (3) 

ai + a 2 

expressed by the lengths of the two semi-axes in the shear plane. For initially spherical 
capsules the deformation parameter D grows with increasing shear rate, while D reaches 
a plateau value for capsules with incompressible membrane area at large shear rate. 

Two main modes of motion of capsules in shear flow can be identified as illustrated 
in figure [2| (i) The capsule shape reaches a stationary state in the laboratory frame. 




Figure 1. Ellipsoidal capsule with semi-axes ai, 02, 0,3 in linear shear flow v°°. The 
shear plane defines the xy plane with shear flow pointing into the x direction. The 
coordinate axes x\, X2, x% point along the principal directions of the capsule where 
the X\ and X2 axes, with a\ > 0,2, are chosen to lie in the shear plane and are rotated 
through the inclination angle \& with respect to the x and y axes. The phase angle 
$ measures the angle between the direction of the shear flow and the direction of a 
tracer particle which initially was on the endpoint of the long semi-axis 01. The phase 
angle </> = $ — * measures the tank-treading motion of the membrane by the angle 
between the position of the tracer particle with respect to its undeformed position. 
The membrane encloses a fluid of viscosity rji and volume V. The fluid outside the 
capsule has viscosity r] . 




Figure 2. Sequences of swinging (top) and tumbling (bottom) of an ellipsoidal 
capsule. The black dot and arrow denote the position and motion of a material 
membrane element. In the swinging motion the inclination angle oscillates around 
a finite positive inclination angle if? while the membrane performs a tank-treading 
motion with monotonously decreasing phase angle <j>. In the tumbling motion the 
inclination angle if? steadily decreases leading to a flipping motion of the capsule and 
an oscillating phase angle <f>. 
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The membrane elements rotate along the fixed shape like the chains of a tank. The 
inclination angle \1/ is then constant, while the phase angle <p grows monotonously. Such 
a motion is called tank-treading in its pure form. A variant of this motion with slightly 
oscillating inclination angle is called swinging and will become important later in this 
paper, (ii) The capsule shape rotates in the laboratory frame as a rigid body would do. 
In this case the inclination angle \l/ grows monotonously, while the phase angle <p stays 
approximately constant. This motion is called tumbling. 

This paper is organized as follows. We give an overview over the theoretical 
methods in section |2j We then summarize the experimental developments in section [3j 
The formation of short-wavelength wrinkles on the surface of polymerized capsules is 
discussed in section |4} We then focus on the dynamics of non-spherical capsules in 
section |5j where we find the predictions of a simple analytical theory at odds with 
numerical evidence. To solve this discrepancy in the main original contribution of this 
paper, we systematically expand the equations of motion in the limit of small capsule 
deformations in section [6j We find that in certain dynamic regimes shape changes, which 
are not covered by the simple theory, play an important role in the capsule dynamics. 
We finish by discussing a few perspectives in section [7} Details of the analytical methods 
and longer derivations are deferred to appendices A-C. 



2. Theoretical approaches 

From a theoretical point of view, the full equations of motion are intractable analytically. 
However, much insight into the physics can be obtained by either studying certain 
asymptotic limits, where the dynamics can be simplified, or by restricting the number 
of degrees of freedom to as few as possible. An example of the first strategy is the 
quasi-spherical expansion of the equations of motion for initially spherical capsules as 



shown in the following section 2.1 The Keller-Skalak model [11 J (described in detail in 
section 5.1 ) follows the second strategy, by artificially fixing the shape of the capsule and 
only allowing tank-treading and tumbling degrees of freedom. In its original form this 
model is only applicable for fluid vesicles although it was intended to describe RBCs. 
It is missing an important energetic contribution to the dynamics, the so-called shape 



memory [12] . The origin of this term is discussed in section 2.2 An extension of the 



Keller-Skalak model including the shape memory effect [T3l ITi] is described in detail in 



section Beyond these two approaches, numerical methods are inevitable. In section 
2.3[ we give an overview of the different numerical approaches that have been developed 
and successfully used to determine the dynamics of vesicles, capsules and RBCs. 



2.1. Initially spherical capsules 

Often the shape of the capsule deviates only little from a sphere. Similarly, in weak 
flow, initially spherical capsules only deviate slightly from their equilibrium shape. 
The highly non-linear equations of motion can then be expanded in terms of the small 
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deformation vector field. One can then take further advantage of the high symmetry 
of a sphere by changing to a spherical harmonic representation. A comprehensive 
treatment is illustrated in |15j . where the displacement field is expanded in vector 
spherical harmonics, which are defined for completeness in Appendix A| since we use 
them in our approach as well. 



In early work, Barthes-Biesel [16J considered the steady state of an initially spherical 
elastic capsule using a small deformation expansion. The membrane is assumed to be 
homogeneous, isotropic, area-incompressible, and shows no resistance against bending. 
To linear order in the deformation only tank-treading motions are possible in the long 
time limit. Only at higher order a transition to tumbling can occur, when the shear forces 
have deformed the capsule sufficiently In the steady state, the capsule is orientated with 
inclination angle — 7r/4 while the membrane performs a tank-treading motion. With 
increasing shear rate, the deformation first increases linearly with quadratic corrections. 
In |17j . the relaxational dynamics of a capsule with arbitrary elastic response was 
obtained still excluding bending energy. There, results have also been obtained in limit 
cases, like droplets (showing only surface tension) or red blood cells (having a local area 
incompressibility). More recent theoretical efforts are reviewed by D. Barthes-Biesel and 
C. Pozrikidis in Chapter 1 and 2, respectively, of [lj. Bending energy was first included 
for capsules by Rochal et al |15j . 



2.2. Shape memory effect: non-spherical capsules 

For a capsule with a spherical equilibrium shape all membrane points are equivalent from 
a mechanical point of view. This is no longer true for an initially non-spherical capsule 
such as an RBC, since the long and short axes are distinguished points on the membrane. 
For a discussion of the consequences for the energetics of a tank-treading motion let us 
assume that the capsule membrane is unstressed in the equilibrium state. As with tank- 
treading motion, we now rotate the membrane by the phase angle <p. Since the local 
geometry of the capsule of each membrane element is different from that of the starting 
position, the capsule membrane will be sheared, requiring a positive elastic energy. 
Consequently, a force aims to restore the unstressed equilibrium state. Only after a 
rotation by <fi = ir the initial and final position of each membrane point are equivalent 
to each other due to the mirror symmetry of the ellipsoidal capsule. The membrane is 
then again unstressed. Any membrane element therefore energetically prefers its initial 
position (or one of the equivalent positions by symmetry of the reference shape). This 
property is called shape memory effect |12j . and plays an important role also in the 
dynamics of red blood cells [HI ESI H2 H3] • 

A direct consequence of the shape memory effect is that capsules perform tumbling 
motions for low shear rates. Here, the hydrodynamic forces are not sufficient to overcome 
the elastic barrier which is present for a hypothetical tank-treading motion. Only for 
higher shear rates the elastic barrier can be overcome and tank-treading motion becomes 
possible. Since the force balance of a capsule is also slightly changed by the shape 
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memory effect, one now observes stable oscillations of the inclination angle \1/ around a 
finite value instead of a stationary tank-treading state. This type of motion has therefore 
been called swinging motion [TJ]. This term should not be confused with the swinging 
mode of fluid vesicles, which is also called vacillating-breathing or trembling [201 |2TI [22] . 

2.3. Full numerical treatment 

For quantitative results beyond asymptotic limits one has to resort to numerical methods 
for capsules in shear flow. Boundary integral methods (BIM) have been used as first 
approaches in handling capsules in flow |231I211I251I2SII2ZIESIESII3DIEIII321I33]. As 
a mixture of global spectral methods [34J and boundary integral methods, the spectral 
collocation [35J and the spectral boundary methods [HS1 EH EH] have been developed 
which use higher order basis functions on a triangulated surface. 

The immersed boundary method [39| SO] can be implemented in any hydrodynamic 
solver. Here, the force density of the boundary is mapped onto a Cartesian mesh. It 
has been combined with solvers of the Navier-Stokes equation |4T| |4"2] and with lattice 
Boltzmann methods [SJ SU SSI SSI El SHI SH]- The front tracking technique |SD] 
can handle two-phase flows and was used with a finite volume method to treat two- 
dimensional capsules [5T] . 

The oscillations of the shape and orientation of non-spherical capsules have first 
been observed by Ramanujan et al [28] who used a boundary spectral method. For high 
and low deformations this method was plagued by numerical instabilities due to the 
degradation of the grid. Further improvement of the boundary element method allows 
the stable simulation of swinging and tumbling of highly flattened capsulses only by 
numerically smoothing the surface [30]. More recent approaches, such as the spectral 
boundary algorithms |3"fi"| [37J EH] use higher order basis functions on the triangulated 
surface to avoid numerical instabilities. Pure spectral methods [34] employ global basis 
functions to parameterize the capsule membrane. The fluid dynamics part is also 
simulated using lattice Boltzmann methods [52, SSI SHI EH]- Recently, front-tracking 
methods have been employed to study the motion of two-dimensional |51] or three- 
dimensional capsules [121 SH] • 

Existing solvers for the dynamics of vesicles, which treat the flow at a continuum 
level, either employ the boundary integral method on a discrete triangulation scheme 
[SU E3J [55] or a phase field or advected field model [56] ETJ [58]. An alternative route 
is provided by stochastic rotation dynamics (SRD) or multi-particle collison dynamics 
(MPCD), where the flow is modeled using effective fluid particles interacting with a 
dynamically triangulated membrane [591 ESI EU EU [62] . 

3. Experiments 

Early experiments |63| EH ESI ES] confirmed the quasi-spherical predictions for the 
deformation parameter D (for not too high shear rates) as well as the decreasing 
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behaviour of the inclination angle ^ for increasing shear rate. Although the inclination 
angle \l/ deviated from the theoretical prediction it /4 even in the limit of vanishing shear 
rates 7 = [64, 66J, experiments together with the quasi-spherical results were used to 
infer values for the elastic parameters of the membrane. 

In contrast to the theoretical predictions for spherical shapes, capsules in such 
experiments [64^, 166] showed oscillations of both the deformation D and the inclination 
angle These oscillations caused by the fact that real capsules deviate from a perfectly 
spherical shape are a direct consequence of the shape memory effect discussed above in 
section 12.21 

In order to observe tumbling motion for experimentally accessible shear rates, 
capsules which show a large deviation form the spherical shape have to be used. Red 
blood cells are an ideal and convenient choice. In early experiments, red blood cells have 
only be examined either in a tumbling state in the blood plasm or in a swinging state 
for high shear rates showing a stationary inclination angle [671 168| |69| [70] . The motion 
for low shear rates and close to the swinging transition was experimentally examined 
by Abkarian et al [14J. Within the swinging regime the frequency of the oscillation of 
the inclination angle is exactly twice the frequency corresponding to the tank-treading 
motion. While a membrane point is advected a full rotation along the membrane, 
the inclination angle is performing two full swinging motions. This relation highlights 
the direct connection between the oscillations of the shape and the motion of the 
membrane. For a decreasing shear rate within the swinging regime, the frequency of the 
oscillation decreases, while the oscillation amplitute of inclination angle increases. The 
numerical analysis of the full equations of motion using spectral methods [34J discussed in 
section 2.3 yields qualitative and quantitative agreement with the experimental results. 
Interestingly, a hysteresis effect at the transition between tumbling and swinging occurs 
when alternatingly increasing and decreasing the shear rate [14J. Here, a single tumbling 
motion seems to mix with multiple swinging motions. 



4. Wrinkling 

Rather than deforming to a slightly ellipsoidal shape, certain polymeric capsules show 
a wrinkling instability while tank-treading in a linear shear flow [UE], see figure [3j The 
wrinkles are caused by a compression-force acting onto the membrane. Although the 
membrane moves along the stationary shape while tank-treading, both position and 
direction of the wrinkels remain constant as they are caused by the compression part 
of the elongational part of the flow. Thin membranes are unstable with respect to 
compressive forces. In a flat geometry, compression leads to Euler buckling on the 
largest possible length scale. To produce small length scale wrinkles on a flat membrane 
one needs a combination of tensile stress along the wrinkles and a non-local geometric 
constraint [71] . 

In the small deformation regime this wrinkling instability was analytically explained 
in [72J using a separation of length scales, since the capsule radius is much larger than 
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Figure 3. Wrinkling of a polysiloxane capsule (~ 343/im) in shear flow as indicated 
by the arrows (adapted from [66 ). 




Figure 4. Wrinkling instability region for 7 = 11. 5s~ 1 (left) and 7 = 17.0s -1 (right), 
respectively. The elastic moduli of the membrane are K = 0.2Nm _1 , k = 1- 10~ 17 Nm. 
In the black coloured region the membrane is still stable with respect to wrinkle 
formations. Brighter colons indicate a larger growth rate of the wrinkles. Initially, 
wrinkling occurs only in a small region around the circumference corresponding to the 
plane of maximal compressive stress. At higher shear rates, this region grows, as does 
the growth rate of each wrinkling mode. The direction of the wrinkles along the lines 
of maximal tension arc indicated schematically (adapted from [72 ). 



the typical wavelength of the wrinkles. Even though the overall dynamics of a capsule in 
viscous flow is inherently non-Hamiltonian, energy considerations dominate the physics 
at the small length scale of the wrinkles. The energy required to form wrinkles of 
amplitude 5r and local wave number k is 

U el ~ Hq + J dAK(5r) 2 /R 2 + a (5r) 2 k 2 + K(5r) 2 k 4 . (4) 

Here, K is the area incompressibility and ctq < the compression produced by the 
hydrodynamic force. This term drives the instability. Large wave numbers k are 
penalized energetically by the bending energy ~ nk 4 . In contrast to a flat membrane, 
buckling on any length scale on a curved membrane costs elastic energy ~ K/R 2 , since 
the area of the membrane is increased (decreased) when a membrane element is moved 



Micro- Capsules in Shear Flow 



9 



outwards (inwards). The resulting tensile (compressive) stress suppresses buckling on 
too large scales. Only when the compressive force is large enough formation of wrinkles 
can reduce the elastic energy in the system. This happens at intermediate wave numbers, 
which can be calculated by minimizing above energy. 

For polymerized capsules, the analytical expression for the critical shear rate above 
which the wrinkling starts reads 

% = < 5) 

and the critical wavelength of the wrinkles becomes 

At 7c, the wrinkles start to form along a circumference which forms an angle of 7r/4 
with respect to the shear direction. At higher shear rates 7 > j c , wrinkling extends over 
larger regions on the capsule, as shown in figure |4j The analytical expression for the 
critical shear rate would allow for an experimental determination of the bending rigidity 
of the membrane which is difficult to obtain otherwise. So far, experimental work seems 
not to have taken up this suggestion. 

5. Reduced model with two degrees of freedom 

We now turn to capsules which are initially non-spherical. 

5.1. Keller- Skalak model 

In the theory of Keller and Skalak [TT], a red blood cell (RBC) is treated within a 
simplified model. Since as an essential ingredient the shape memory effect of RBCs is 
missing within this simplified description, the model does not capture essential parts of 
the physics. Yet, for vesicles the Keller-Skalak model delivers quite a good description. 
Since it is also the basis for the further modelling of capsules, we briefly recall its main 
features. The three-dimensional RBC or vesicle is assumed to have a fixed ellipsoidal 
shape 

(xi/ai) 2 + (x 2 /a 2 ) 2 + (x 3 /a 3 ) 2 = 1, (7) 

see figure [T] 

The (Cartesian) components of the undisturbed shear flow are (7?/, 0,0). In order 
to mimic a tank-treading motion of the membrane, the velocity field at the membrane 
is assumed to have components 

v = w*f (-(oi/a 2 )x 2 , (02/01)0:!, 0) , (8) 

where u^ s is a parameter having the dimensions of a frequency. The energy supplied 
by the external flow has to be balanced with the energy dissipated inside the RBC. 
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* = - ( i + -i^p-j + cosm . (9) 



Likewise the torques acting on the RBC have to balance. The equation of motion of the 
RBC derived from this torque balance reads |TTJ 

2 a? + a| / af + a|2 
The dissipated energy must come from the work performed by the external flow, which 
leads to the equation 

nJitf = Vo(f2<p 2 + / 3 70cos(2*)) , (10) 
where the geometric quantities depend only upon the ratios of the axis lengths dj and 
are defined by 

1 fax a 2 



cti = di (a^Os) 1 ^ 3 for i = 1, 2, 3 , z\ 

2 \a 2 cl\ 



Jo 

Both balance laws can be simplified to 

4l = ^(A + Bcos(2^)), (12) 
= 7 Ccos(2^) (13) 

with 

A =-1/2, (14) 

d _ Q1Q2/3 + (/2 - fiVi/Vo)(a 2 i ~ of) nt -N 
" 2(a 2 + a|)(/ 2 -M/r /o ) ' 1 ° j 

g = ~ f t / • ( 16 ) 

J2 - flVi/Vo 

The phase diagram for the dynamics of RBCs or vesicles is rather simple in the 
Keller-Skalak model. The transition between tumbling and tank-treading is independent 
of the shear rate. For given half axes a,, there is a critical viscosity contrast (r)i/r) ) c 
given by \B\ = \A\ = 1/2, below which the RBC or vesicle performs a tank-treading 
motion. For larger viscosity contrast, it undergoes a tumbling motion. The transition 
from tumbling to tank-treading is thus independent of the shear rate 7 within the 
Keller-Skalak model. 

For low viscosity contrast i]i/r] < (r)i/r) ) c , i.e. \B\ > 1/2, one obtains a steady 
tank-treading angle 

v^arccosQ (17) 

and a steady tank-treading frequency 

^ = i^- (18) 

For high viscosity contrast r]i/r] > (r)i/r) ) c , i.e. \B\ < 1/2, one obtains a steady tumbling 
motion with mean tumbling rate 

(*> = Hm 1 [ T ^(t')dt' = -ly/T=4B*. (19) 
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5.2. Skotheim-Secomb model 

The Keller-Skalak model in its original form can only correctly describe the motion of 
vesicles, as it does not take into account the shape memory effect. Abkarian et al jH] 
and Skotheim et al [13J extended the Keller-Skalak model to remedy this deficiency. 

The membrane with its elastic property serves as an energy storage and induces an 
additional torque onto the capsule. The elastic energy E((f>) required for a hypothetical 
tank-treading motion must be positive and periodic in the phase angle <p with period 



7r, as argued in section 2.2 In [Hf [13], the energy was modeled as a simple harmonic 
function 

E{<f>) = £ o sin 2 0- (20) 

A Keller-Skalak derivation balancing energy and torque leads to two equations of motion 
for the two degrees of freedom, the inclination angle and the phase angle. 

4f(t)=A(t)+Bcos(2V(t)), (21) 
0(t) = C (cos (2*(f)) - sin (20(f)) J (22) 



with constants 



.4(7) ^ ( . o(/)V (23) 




B s 25J+3)* (24) 

C = - f . 3 7- (25) 

f 2 - hmho 

In dimensionless quantities of shear rate x, viscosity ratio A and eccentricity a 

X = ^ 7, (26) 

+ ^ (27) 

a 2 - a 2 

a = arctan — -, (28) 

2ai<22 

the equations of motion can be cast in simple form [73] 

4>(r) = - cosa0(r) - A(l - sin a cos2^(r)) , (29) 
0(r) = - (x~V)sin2(/>(T) + cos2^(r)). (30) 

Here, where r is a dimensionless time defined differentially as 

27 

(It = -^dt. (31) 
A 

Solving the equations of motion within the reduced model for arbitrary parameters 
has been performed numerically [H f [T31 173] . Kessler et al [73] have obtained analytical 
solutions in the quasi-spherical limit and, using of the method of asymptotic expansion 
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Figure 5. Grey scale plots of the normalized mean tumbling rate within 
the reduced model for two different ellipsoidal reference shapes with eccentricity a, 
a = corresponding to a spherical shape within the shear plane. Dark grey colour 
corresponds to a vanishing mean tumbling rate uj tu — 0. while light grey corresponds 
to pure tumbling u tu = 1. 



[75| . both the stable trajectories and the slope of the transition between tumbling and 
swinging. 

For small shear rates 7 the capsule always tumbles as the hydrodynamic forces are 
too weak to overcome the elastic barrier. For increasing shear rate the transition to the 
intermittent motion takes place. By further increasing the shear rate, the transition to 
swinging occurs for small and moderate viscosity contrast A. 

The phase diagram (see figure [5]) of this reduced model shows the regimes of the 
three possible motions. Either the capsule tumbles, swings or performs an intermittent 
motion with mixed tumbling and swinging periods. The normalized mean tumbling rate 

<*> s lim fo^m/T 



(32) 

(*> + <0) T^ff*(t)dt/T + ff4>(t)dt/T 

serves as an order parameter to distinguish between the possible motions [73]. Pure 
tumbling motion is characterized by u tn = 1, pure swinging by u tn = 0, while 
for intermittent motion the mean tumbling rate u tu takes values between and 1. 
Noguchi has pointed out that u tn reaches several plateau values between and 1 in the 
intermittent regime |74| . 



5.3. Comparison with spectral method 

To check the validity of the reduced model, its predictions have to be compared with the 
full non-linear dynamical motion obtained by numerical methods. The first systematic 
numerical phase diagram for a large range of shear rate and viscosity contrast was 



constructed using the spectral method described in section 2.3 [34J. Using this method 
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Figure 6. Comparison of phase diagrams deduced by the spectral method [M] 
(symbols and full line) with the one following from the reduced model (dashed lines). 
Phase diagrams are functions of the viscosity contrast r\ jr\i and the inverse shear rate 
(Rf]oi I lA^ 1 1 where /j, is the shear modulus of the capsule. For low shear rates the 
capsule tumbles while for higher shear rates the capsule performs swinging motions 
for a not too high viscosity contrast. In between these two regimes, there is either a 
transient regime or intermittent regime deduced by the spectral method or the reduced 
model, respectively. Adapted from |34j . 

tumbling motions for capsules slightly deviating from the spherical shape were observed 
numerically for the first time. 

The phase diagram obtained from the spectral method as shown in figure [6] shows 
qualitative agreement with the predictions of the reduced model. The reduced model 
correctly predicts a swinging and tumbling motion and a transition for increasing 
shear rate. In particular for small viscosity ratios, the slopes of the tumbling-swinging 
transition agree quantitatively, while the absolute position of the transition differs by 
about 20 percent. 

Yet, in the intermittent regime of the reduced model, only transient motion was 
found in the full model as can be seen in figure [7j Here, the capsule initially tumbles 
with an steadily decreasing inclination angle while the phase angle is oscillating. 
After a few tumbling motions, the capsule starts to perform stable swinging motions, 
where the inclination angle \I/ oscillates around a stable mean value, while the phase 
angle <fi is steadily decreasing. Independent studies using two different solvers confirmed 
this finding: A transient motion relaxing from a tumbling motion into a stable swinging 
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Figure 7. Transient motion from tumbling to stable swinging motion. In this 
parametric plot the transient trajectory is shown as a function of the inclination angle 
*S> and the phase angle <f> where the arrow denotes the direction of time. Initially, the 
capsule tumbles with oscillating </> and steadily decreasing , but finally ends up in 
a swinging motion, where <!' begins to oscillate around a stable value. Adapted from 



motion was always observed in the intermittent regime |4T1 142] . We therefore prefer to 
call the disputed regime transient rather than intermittent |34| . 

Close to the transition to tumbling, the oscillation of the inclination angle \1/ is very 
strong, surprisingly leading to temporarily negative values of the inclination angle, as 
first noted by Kessler et al j34]. Here, the shape changes are very pronounced. This 
effect was somewhat unfortunately named breathing by Bagchi et al |42]. It should not 
be confused with the vacillating-breathing motion for vesicles |20| because it does not 
show pronounced deformations perpendicular to the shear plane. Moreover, the mean 
inclination angle is not close to zero but finite. Since in all studies shape changes are 
highly pronounced in the transient regime (see (42] in particular), this reduced model 
has to be expanded to allow for shape changes to unambiguously clarify the status of 
the disputed intermittency. 

6. Systematic quasi-spherical expansion 

From the numerical findings it is clear that allowing shape changes is one important 
ingredient needed to solve the open intermittency question. Noguchi [74] incorporated 
shape changes within the reduced model heuristically, following extensions of the Keller- 
Skalak model for vesicles [21] . Besides inclination and phase angle, a shape parameter is 
used as a third degree of freedom. The elastic energy due to shape changes is numerically 
calculated by stretching a biconcave membrane (with a Skalak energy law). The equation 
of motion of the shape parameter is deduced by a deformation from the sphere, while 
the equations of motions for the two angles are obtained by a Keller-Skalak derivation 
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[llj. The overall phase diagram of the extended model was qualitatively no different 
from that of the reduced model. In particular, the intermittent motions survived. A 
temporarily negative inclination angle \1/ is also found in the expanded model. The effect 
of the bending energy was studied using the same heuristic approach |76| . 

6.1. Equations of motion 

We would like to replace this ad-hoc treatment of shape changes by a systematic 
expansion of the full equations of motion. Specifically, we expand the equations of 
motion of the capsule for nearly spherical shapes. While the deformation of the capsule 
shape from a sphere stays small, the displacement of a membrane element during a 
tank-treading motion is usually large and cannot be treated as a perturbation. The 
derivation of the elastic energy and the equations of motion is straightforward, but 
would distract from the main physics. Details are therefore presented in appendix B 
and appendix C. From the elastic energy, we derive the elastic forces, which must be 
balanced by hydrodynamic forces. The hydrodynamic forces have contributions from 
the externally imposed flow and from the induced velocity field, which must be solved 
dynamically. No-slip boundary conditions at the membrane and at the external walls 
uniquely determine the dynamics. 

Rather than specializing to linear shear flow, we consider general linear external 

flow 



Here, s is a measure of the strength of the elongational flow component and 2lo is 
the vorticity. Linear shear flow of the form v°° = A fye x corresponds to s = j/2 and 
uj = —7/2. In a microfluidic four-mill device, which was used to study vesicle dynamics 
quantitatively |77j . it is possible to control u and s independently. 

For motions with the same mirror symmetry with respect to the xy-plane as the 
external flow we derive the equations of motion for the inclination angle phase angle 
0, and a shape parameter (3. The shape parameter j3 G [0, 7r/2] is a measure how much 
of the excess area is stored in the m = 2 modes. For (3 = the capsule is a prolate 
rotationally symmetric ellipsoid with the long axis in the vorticity (z) direction. For 
(3 = 7r/3 the capsule is again prolate rotationally symmetric, but with the long axis in 
the shear (xy) plane. (3 is connected to the Taylor deformation parameter defined in ^ 
to lowest order via 



The corresponding shape parameter of the equilibrium shape is denoted with a hat f3 
(without loss of generality ^ = $ = 0). The equations of motion as derived in the 
appendices read 




(33) 




(34) 




A 



(35) 



S' 1 sin 20 -cos 2$ 



(36) 
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$ = - S~ 1 cot/3sin/3 + cos/3(^ 1 cos20 + sin2#), (37) 



with 



/30tt 



T 



A 32 + 23^/77,, 



(38) 



S^J^^Ks, (39) 
V A //sin/3 V ; 

A=-^A(32 + 23W<. (40) 

Here r denotes a dimensionless time and the dot the corresponding derivative. S is a 
measure of the strength of the elongational flow, and A measures the relative strength 
of the rotational flow. In linear shear flow (s = — u> = j/2), A depends only on the 
viscosity contrast and material capsule parameters, while S is proportional to the shear 
rate 7. 

6.2. Numerical results 

Before we proceed to analyze these equations in more detail, we summarize the 
findings based on numerically integrating these equations for an initially axisymmetric 
equilibrium capsule with /3(0) = ft — 7r/3. In figure pi the mean shape parameter (/3) in 
the long time limit is plotted colour coded as a function of the viscosity contrast A and 
the elongational shear rate S. In this "dynamical phase diagram", one can distinguish 
three regions: (i) A swinging region, roughly bounded by A < 1 and S > 1. The mean 
tumbling rate u tn ~ vanishes in this region, while the shape parameter /3 ~ 7r/2 is 
maximal. A typical trajectory over time is shown in figure [9] (ii) A tumbling region 
roughly bounded by A > 0.6275. Here, u tn ~ 1 takes the highest possible value, while 
the shape parameter oscillates around the equilibrium one, f3 ~ /3, in the developed 



tumbling regime. A typical trajectory is shown in figure 10 (iii) A transient region in 



between. Here, the shape parameter (3 ~ is small, while the tumbling rate m tu ~ 



vanishes in the long time limit. A typical trajectory shown in figure [TT] reveals that the 
capsule initially tumbles for a while, before it assumes a stationary swinging state. The 
initial tumbling time increases with increasing A. 

The shape parameter can therefore act as an order parameter for distinguishing 
the dynamical phases. A cut through the dynamical phase diagram at S — 10, initially 
starting at the equilibrium shape (3(t = 0) = (3 is shown in figure 12 We see that the 



swinging to tumbling transition at A ~ 1 is continuous, while /3 jumps at the transient 
to tumbling transition at A ~ 6.27, like in a first order phase transition. 

6.3. Analytical results 

After having presented the numerical phenomenology, we analyze the equations of 
motion analytically and make the connection to the reduced model of section [5j where 



the shape (3 of the capsule is held fixed. If we ignore the dynamics (37) of the shape 
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Figure 8. Dynamical phase diagram obtained by numerically integrating ( 35 37 1 
for an initially axisymmetric equilibrium capsule (/3(t = 0) = j3 = 7r/3). The 
mean shape parameter (f3) is shown in a grey scale plot as a function of the shear 
rate S and the viscosity contrast A. Dark grey corresponds to the minimal value 
(/?) = 0, while light grey denotes a maximal value (j3) = ir/2. Three regimes can be 
distinguished: swinging, tumbling and a transient motion towards swinging. The small 
circles correspond to figures [&ffl"H while the cut at S = 10 through the phase diagram 
is shown in figure [12] The vertical dashed line corresponds to the analytically obtained 
phase boundary at A = 1 between swinging and transient motion in the limit S — > oo 
with constant A. The inclined dashed white line corresponds to the phase boundary 
between tumbling and transient motion in the limit S, A — > oo with constant ratio 
A/5 ~ 0.627. 




tA/tv tA/tt 



Figure 9. Shape parameter f3 (left), inclination angle ^, and phase angle </> (right) as 
a function of dimensionless time tA/t: within the swinging regime (/3 = 7r/3, A = 0.5 
and S — 6.0, see figure [8]). The shape parameter j3 oscillates around a value (f3) close 
to 7r/2, implying an almost maximally extended shape within the shear plane. The 
inclination angle ^ oscillates around a constant positive value (\&), while the phase 
angle (f> linearly decreases with superimposed oscillations (not seen in the plot due to 
scale) . 
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Figure 10. Shape parameter f3 (left), inclination angle "J, and phase angle <f> and as 
a function of dimensionless time tA/tt within the tumbling regime (f3 = tt/3, A = 5.0 
and S — 6.0, see figure |8j. The shape parameter /? oscillates around a value close to its 
reference value j3 — 7r/3. The inclination angle ^ linearly decreases with superimposed 
oscillations, while the phase angle tfi oscillates around a constant positive value ((/>). 
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Figure 11. Shape parameter f3 (left), inclination angle fy, and phase angle <fi (right) as 
a function of dimensionless time tA/t: within the transient regime (f3 — 7r/3, A = 2.5 
and S = 6.0, see figure[8]). Initially, tumbling motions occur where the inclination angle 
^ linearly decreases with superimposed oscillations, while the phase angle (j> oscillates 
around a constant positive value. On a longer time scale, the shape parameter f3 relaxes 
to a smaller value, where the capsule is close to circular in the shear plane and strongly 
extended in the vorticity direction. After the relaxation, the motion changes into a 
stable swinging, where oscillations are superimposed to a constant mean inclination 
angle (fy) and a linear decreasing phase angle (j>. 



parameter and compare the remaining equations (35, 36) with the corresponding ones 
in the reduced model, (29, 30), we see that up to a rescaling in time both equations are 
identical in the quasi-spherical limit (a — > 0). We identify 

rsK^TQs/sin^), (41) 

A =Asin(/3), (42) 

X =S. (43) 

Thus, if the dynamics of the shape parameter is slow compared to that of the angles, 
we can read off the fast capsule dynamics in the phase diagram of the reduced model, 
where the effective viscosity ratio A changes slowly in time. We will see that for strong 
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Figure 12. Cut through the phase diagram of figure [8] The numerically obtained 
mean shape parameter (f3) is plotted as a function of the viscosity contrast A at 
constant shear rate S — 10 for an axisymmetric capsule starting out of its reference 
shape /3(0) = f3 = ir/3 (black line). The grey lines correspond to asymptotic analytical 
results. For low A < 1 the capsule swings with maximally extended shape within the 
shear plane. For high values A > 6.27 the capsule performs tumbling motions with 
shapes close to the reference shape. In between there is the transient regime where the 
extension of the shape in the vorticity direction strongly increases with increasing A. 
The arrows denote the flow towards the long-time stationary value. 



flows the shape changes indeed on a much slower time scale than the angles. This 
separation of time scales is the key to understanding the dynamic phase behaviour of 
capsules in the quasi-spherical description. 

6.3.1. Swinging /transient transition, S ^> 1, A fixed: We start our analysis by taking 
the limit of strong flow S — > oo, while keeping the relative strength of the rotational 
flow A constant. The leading order equations then become 

^ + = - A, (44) 

(9 + A)sin/3 = cos 2^, (45) 
$ = cos (5 sin 2tf. (46) 

We can identify two dynamical regimes. 

Swinging: (A < 1): We now look for stationary tank-treading solutions for the 
inclination angle \I/ = and the shape parameter (3 = /3q. Assuming that sin2\l / 
stays positive (0 < < 7r /2), the shape parameter approaches the value f3 = n/2 
asymptotically. The remaining two equations then reduce to the Keller-Skalak form 
^ = —A + cos2 , I> with the stationary solution for A < 1 

^ w = ^arccosA. (47) 
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The next order in the inverse shear rate yields oscillating perturbations of the inclination 
angle 

i£ sw = I arccos A + cos2Ar + 0(1/ 5 2 ), (48) 
while the shape parameter is shifted with second order oscillations 

r=-~ m ^ + f $ sin2Ar + 0(l/5 3 ). (49) 



These results correspond to the numerically observed swinging dynamics. In figure 13 



the numerically calculated long term mean inclination angle (ty) is compared to the 
corresponding first term on the right hand side of (48 ) as a function of A. The amplitude 
of the oscillatory component, the second term on the right hand side of (48 ), is compared 
to numerically obtained data in figure 14 The amplitude of the second order oscillations 



of the shape parameter (3 finally is shown in figure 15 



Transient (A > 1): For A > 1 (44-46) admit a different stationary solution given 

by 

#Jf = and ffi = arcsin(l/A). (50) 

These values define the stationary state in the transient regime, which is only reached 
after initial tumbling motions. The first order calculations add oscillatory terms 
according to 

tr = -Ar TOt/3 - lf~ 1 - cos2Ar + 0(1/ S 2 ), (51) 

25 V / A T ^I 25(A 2 + 1) v 1 h y ' 

tr = c ot/3 3A 2 - 1 0(l/5 2 ), (52) 

2SVA^1 2S(A 2 + 1) v 1 h y 1 

f3 tr = arcsinl/A + sin2Ar + 0(1/ S 2 ). (53) 

Again, the mean value of the inclination angle and the oscillation amplitudes of the 



inclination angle and shape parameter are compared to numerical data in figures [13 - 15 



6.3.2. Transient /tumbling transition, S, A 3> 1,A/S fixed: In order to discuss the 
tumbling regime, we must consider a different asymptotic limit. The numerical results 
indicate that the transient to tumbling transition occurs for large shear rates at 
fixed ratio A/5. In this limit, two timescales must be considered. On average the 
parameters vary on a slow timescale r, but the fast oscillating trigonometric terms 
produce oscillations on the much faster timescale Ar. We seek an expansion of \1/ ((f), (3) 
into powers of the small quantity 1/A, ^ = + ^i + • • •• The time derivative of the 
small oscillating corrections are multiplied by A and can thus contribute to the order 
one on the left hand side of the equations of motion. The systematic method to analyze 
processes with different time scales in a perturbation scheme is described in |75| . 
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Figure 13. Mean inclination angle (if?) as a function of the viscosity ratio A for 
j3 = 7r/3 and S = 10. The black line corresponds to the numerically obtained values, 
while the grey lines show analytical results. The mean inclination angle decreases from 
7r/4 (corresponding to the direction of the elongational part of the shear flow) at A = 
in the swinging regime. At the transition to the transient regime, the mean inclination 
angle is close to 0. In the transient regime, (if?) asymptotically reaches 0. 

0.05 1 — , — , — , — , — i — , — , — , — , — i — , — , — , — , — i — , — , — , — , — i — , — , — , — , — i — , — , — , — , — | 




Figure 14. Amplitude A^f* of the oscillations of the inclination angle $ as a function 
of the viscosity ratio A for f3 = tt/3 and S = 10. The black line corresponds to the 
numerically obtained values, while the grey lines show analytical results. 



Long time swinging behaviour: First, we seek to rediscover the stationary tank- 
treading solution in the transient regime. We expect ~ —A to leading order. Since 
only terms of order one occur on the right hand side of (36), sin/3 must be small of 
order 1/A on the left hand side. Taking this into account, we recover the tank-treading 
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Figure 15. Amplitude A/3 of the oscillations of the shape parameter j3 as a function 
of the viscosity ratio A for f3 = ir/3 and 5 = 10. The black line corresponds to 
the numerically obtained values, while the grey lines show analytical results. In the 
swinging regime, the oscillations show up only in second order in 1/5, while in the 
transient regime, the amplitudes are of order 1/5. Shape changes of the swinging 
capulse are small in the swinging regime and pronounced in the transient regime. 



solution (arcsin(l/A) «1/A) 

#7 = and A 



sw 




1/A. 



(54) 



Expanding the equations of motion ( 35 - 35 ) to the next orders in 1/5* gives the oscillatory 
corrections 

3 



A s 



2S 



cos 2Ar + 0(1/5" 



- Ar- — cos2Ar + (9(l/5 2 ), 
l/A + -^sin2Ar + C(l/ 1 S 3 ). 



(55) 
(56) 
(57) 



Note that the swinging motion is linearly stable for arbitrarily high values of A. 

Tumbling: During tumbling motion the inclination angle # grows very fast, to 
leading order # ~ —A. We assume a slow dynamics of the shape /3q(t), superimposed 
by small but fast oscillations Ai(At). From the leading order equations of (35-37) we 
obtain 

, . 1 



2 A sin A 



sin(2A7 



ft ( r ) = ~^2A~ COS 
The equations at order 1/A give 

• 1 sin 20o / cos 2Ar 

S sin/3o \sinAo + Ai 



(58) 
(59) 

(60) 
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11 + sin 2 /3 
4A sin 2 fa 



S sin/3 

1 1 1 

cot (3 sin (3 + — cos (3q cos 20 o + — (cos (3 sin 2^) , 



(61) 



where (• 



S ' " u S S 
•) denotes averaging over one period of the fast oscillations. The average of 



the cross term (cos/3sin2\I/) on the right hand side of (61) is zero, since the oscillations 
of and \1/ have a relative phase shift of n. 

Stable tumbling regime: Equations (60-61) admit a stationary solution determined 
by the conditions 

1 + sin 2 fa S 



sin 20o 
cos 20 o 



sin/3 4A' 



tan /S 
tan (3 



(62) 
(63) 



The solution (3(A/S) of this system, which is shown in figure 12, has two branches. The 
upper branch was found to be numerically stable. For A/S — > oo the stable branch 
approaches /?o — > /3 (with 0o 0). On the unstable branch we have (3q — > (with 
0o - > 71 /4) as A/S -> oo. We can avoid multi-valued functions by instead writing A/S 
as a function of f3 



A 

S 



(3 - cos 2#; u ) tan0 



tu 



8 sin A tu 



tan 2 /% u 



(64) 



j Q \J tan (3 — tan t j 

The tumbling solution only exists above a critical ratio A/S\ crit ~ 0.627. This finding 
corroborates the first order phase transition found numerically. An analytical expression 
for A/S\ crit can be derived by minimizing the right hand side of (64). However, the 
resulting expression as a function of (3 is too involved to be reproduced here. 

The flow field of the shape parameter (3 is shown in figure 12 for S = 10. Below 
A ~ 6.27, there is only one (oscillating) stationary solution for (3, which is the stable 
transient motion. Above A ~ 6.27, both the stable and unstable tumbling branch and 
the swinging solution coexist. If we start with an initial capsule shape (3(0) above the 
value (3q u (A/S) on the unstable branch, the capsule tumbles in the long time limit. 
Starting with an initial capsule shape /3(0) below the value (3q u (A/S) on the unstable 
branch, the capsule will swing. To check the coexistence of both motions, we solved the 
dynamics now starting with the initial shape /3(0) = arcsin(l/A). With this carefully 
chosen initial condition, transient motion is indeed stable for arbitrarily large values of 
A. 

Initial tumbling in the transient regime: Equations (60-61) can also lead to 
more insight regarding the transient regime. Since no tumbling solution is stable 
for A/S < 0.627, the shape parameter monotonously decreases until it reaches /3 = 
arcsinl/A w 1/A. While (3(t) is of the order one, the slow shape dynamics is given by 
(60 61), which can be easily solved numerically. It is, however, instructive to consider 
the parameters 1 < A < S well inside the transient regime away from the phase 
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boundaries. Here, the equations of motion simplify to 



1 l + sin 2 /3 
J 
1 



(3^= - ^cot/3sin/3 . (66) 

We see that two slow time scales become relevant: The phase angle <fi changes 
monotonously on a typical time scale A, while the shape relaxes on the much longer 
time scale S. During this time the capsule tumbles. The decreasing phase angle 



during the initial tumbling motion is also visible in figure 11 



7. Perspectives 

In this paper we have summarized the current state of research concerning the dynamics 
of capsules in hydrodynamic shear flow. We have reviewed and categorized the 
phenomenology observed in experiments and in numerical simulations. In the review 
part, we focused on the small deformation of initially spherical capsules in shear flow 
and on the formation of wrinkles on polymerized membranes. For initially non-spherical 
capsules, the shape memory effect drastically alters the dynamics, leading to a swinging 
motion. A reduced model with fixed shape explains the swinging behaviour at low 
viscosity contrast, but is at variance with numerical results at large shear rate. We 
have solved this discrepancy by systematically expanding the equations of motion for 
small deformations. Shape changes were found to be paramount to understanding the 
dynamics. The disputed intermittent regime was found to be an artefact of artificially 
fixing the capsule shape. In the relevant parameter range the motion is rather a transient 
one as suggested by numerical evidence. We have obtained analytical expressions for the 
phase boundaries at strong flow and have identified the swinging to transient transition 
as being of second order, the transient to tumbling transition as being of first order. 
In the tumbling regime a swinging motion is also stable for carefully chosen initial 
conditions. 

Until recently, all theoretical and experimental works have focused on soft objects 
in time-constant flow. The dynamics within the reduced model with constant shape was 
completely solved analytically for arbitrary time- dependent flow in the quasi-spherical 
limit [M]. Here, the shear flow oscillates around a finite positive value. With suitable 
modulation of the shear rate, resonant tumbling motion can be excited, even though the 
mean shear rate lies inside the swinging regime for constant flow. First experiments in 
time- varying external flow considered red blood cells [18], where the shear rate oscillates 
with a vanishing mean value. The focus of this study was on the deformation and 
relaxation of the shape. The oscillating shape showed a time lag with respect to the 
shear stress of the flow. The shape changes were characterized by a rapid deformation 
followed by a slow relaxation. Within the expanded reduced model, Noguchi [75] found 
a qualitatively similar behaviour, deducing that in the experiment the red blood cells 
have been in the low-frequency swinging regime. For vanishing mean shear rate, there is 



Micro- Capsules in Shear Flow 



25 



a rich phase behaviour consisting of swinging and tumbling motions and shape changes. 
For finite mean shear rates, the shape changes also lag behind the oscillating shear stress. 
These results imply that shape changes might be important in time-varying flow. 

How the motion in time-dependent shear flow is altered by allowing shape changes 
is an open question worthwhile pursuing. We expect particularly interesting effects in 
the tumbling regime, where two modes of motion are simultaneously stable with their 
own basin of attraction. When changing the flow parameters periodically, hysteresis 
might play an important role. 

Shear flow above a substrate to which a capsule adheres could lead to dynamical 
unbinding, or detachment. For vesicles, this effect has been studied both theoretically 
[79| 180] I8~T] and in experiments [82], [83], [M]. For capsules, the equilibrium shape of such 
a bound configuration has been calculated as a function of adhesion energy and material 
parameters [85, 86J. It would now be interesting to investigate how shear flow affects 
such shapes and how shear-induced detachment of capsules differs from that of vesicles. 

A related, but different open area of research is the influence of thermal fluctuations. 
While thermal noise is energetically not as relevant for capsules as it is for vesicles, 
it is known to influence phase transitions. Thermal fluctuations can be added to 
the equations of motion for capsules as has been done for three-dimensional quasi- 
spherical vesicles [87, SB] [EH] EI] and two-dimensional quasi- circular vesicles [90, 91J. In 
Steinberg's group [92J, thermally induced fluctuations in shape and inclination angle of 
a giant vesicle in a stationary tank-treading regime was investigated as a function of the 
membrane area and viscosity contrast. 

Finally, it looks promising to treat the stationary motion of capsules and vesicles 
as examples of non-equilibrium steady states from the perspective of stochastic 
thermodynamics |93| 194] . One could then ask for the entropy production in such states 
and its behaviour at these dynamical transitions. Likewise, the response of swinging 
and tumbling motion to additional small perturbations could then be expressed by 
correlation functions in the corresponding steady state |95| . Within such a perspective, 
vesicles and capsules could become two-dimensional paradigms for exploring these 
novel concepts which so far have mainly been studied for zero-dimensional and one- 
dimensional objects by using, respectively, colloidal particles and polymers. 

Appendix A. Vector spherical harmonics 

A complete orthonormal basis for three-dimensional vector fields on the sphere is given 
by the vector spherical harmonics |96| which are defined as 






1 




(A.3) 
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where Y l m (9,4>) are the standard scalar spherical harmonics (7 = 0,1,... and m = 
—I, + the operator 

V = e e d e + —^d^ (A.4) 
sma 

is the surface gradient on the sphere, and the scalar product of complex vector fields a 
and b is defined as 

/•27T rn 

(a,b)=/ / a*-bsin#d#d0 (A.5) 
Jo Jo 

with index * denoting complex conjugate. 
Appendix B. Elastic energy 

Each membrane point can be labelled by two Lagrangian coordinates (0r,0r), or 
equivalently by points x(# R , R ) on a hypothetical reference sphere, if (6*r, </>r) are 
interpreted as spherical coordinates. Here, 

x(# R , 0r) = -R[sin(6» R ) cos(0 R ), sui^r) sin(0 R ), cos(6» R )] T (B.l) 

are the Cartesian coordinates of a point on the surface of a sphere with radius R with 
given spherical coordinates. The actual membrane configuration is then constructed in 
two steps: First, the reference sphere is rotated by a three-dimensional rotation <fr = $n 
by an angle <3? around an axis of rotation n. The position of the membrane point is then 
given by D$x(#r, 0r) = x(#,0). Here, D$ represents the Cartesian rotation matrix. 
The analytical expressions for the rotated spherical coordinates (8, <fi) as a function of 
the original coordinates (#r, 0r) and the rotation angle $ are involved. Only in the 
special case where n = e z they simplify to 9 = #r, <p = 0r + $. 

In a second step, the points of the rotated sphere are displaced by a deformation 
vector eRu(8, 0). The final position of a membrane point is 

r(0R,0 R )[$,u] = D*x(#R,0 R ) + ei?u(M). (B.2) 

We have factored the capsule radius out of the definition of u to make u dimensionless. 
The pre-factor e formally a small parameter for our perturbation expansion. Its specific 
value will be chosen below for a convenient normalisation of u. The equilibrium shape 
of a capsule is given by a specific deformation u (we can set $ = without loss of 
generality) . 

The elastic energy of such a deformation is modelled using an Evans-Skalak 
constitutive law with locally and globally incompressible membrane [lj. In the quasi- 
spherical limit the energy can be written as 

H[u,&]= ii J dAe:e + J dAa. (B.3) 

The local tension a is a Lagrange multiplier to ensure local and global area 
incompressibility. The strain tensor 

e(0 R ,0 R ) = l - [g[u](M) " £*g[u](M)] (B.4) 
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measures the difference in the metric tensors g between the equilibrium and the actual 
configuration. Note that the metrics are evaluated at the same material point 0r), 
which are in general at different Cartesian positions in both configurations. One metric 
tensor must therefore be moved to the other position. This is performed using the 
rotation operator D$>, which transforms the tensor components accordingly. We can 
simplify the calculation of the elastic energy considerably by taking advantage of the fact 
that the strain tensor is linear in the deformation in the quasi-spherical approximation. 
The strain tensor and the corresponding elastic energy "Ho[ u ] f° r a deformation from a 
spherical equilibrium shape (u = 0) without area constraint has already be calculated 
in the literature [T5| ITS] . The corresponding expression for a deformation from a non- 
spherical equilibrium shape is then simply %o[ u — D&u]. 

To take advantage of the spherical symmetry of the reference sphere, we switch to 
a vector spherical harmonics representation of the deformation field 



u 



"oo 



47r Im 



Here and in the following, the sum over I and m runs from I > 1 with —l<m<l. 
The vector spherical harmonics are defined in Appendix A Expressions for the volume, 
area, and bending energy in terms of the deformation amplitudes uj m are derived in 
|89j, while expressions for the elastic energy against shearing and stretching in the case 
of spherical capsules can be found in [15J. 

To obey volume conservation, the constant Woo i s a function of the higher mode 
amplitudes and can be eliminated. The excess area is then given by 

A c 2 ST^ i Y + 2 )(^ ~ jj _ g2 PlJ |2 m c\ 

/X - e 2^\ u lm\ 2 = -^2^*l\ U lm\ ■ l B - b ) 

Im Im 

The shear modes uf m completely decouple from the other two deformation modes 
and are omitted in the following. Taking advantage of the local area constraint 
2{ u Ym ~ D&uY m ) = yl(l + l)uf m we can write the elastic energy as 

U[u, <&] = e*R^ E 4( C, + ^ 2 V L - D*uJJ (B.7) 

Im V / 

+ e 2 i? 2 |E*rj 2 . 

Im 

Here, a denotes the homogeneous tension a = J dAa/A ensuring global area 
compressibility. 

Although e is formally a free parameter signifying the order of magnitude of the 
terms, we chose to set 



e = y/A/2. (B.8) 
With this choice the normalization of the U; m reads 

\i:FH m \ 2 = l. (B.9) 

Im 



Micro- Capsules in Shear Flow 



28 



When only / = 2 modes are excited, this simplifies to J2m \ u 2m\ 2 = 1- 

In the spherical harmonic representation the rotation operator reduces to the 
Wigner rotation matrices 

£»«L = E^m(*K»'- (B.10) 

m! 

The elastic force in radial direction is given by 

f Fel _ 9H[U] 

Jim = ~ , ma v* (B- 11 ) 



/ie4(/ 2 + /-2) 
E /(/ + !) 



where the global area constraint contributes a tensile force. Additional forces ensure 
local area incompressibility under the capsule motion. 

Appendix C. Equations of motion 

Since Stokes' equations are linear, any solution can be written as the superposition of 
different flow fields. Moreover, we can chose the boundary conditions at the membrane 
and at infinity of each particular solution to our advantage, as long as the overall solution 
matches the correct boundary conditions given by the force balance and the external 
flow. We can thus separate the flow fields induced by an elongational and rotational 
part and flow induced by the membrane forces and torques 

v(r) = v ind (r)+v s (r)+v w (r). (C.l) 

For the rotational velocity field v w we take the corresponding component of the 
external flow throughout 

v w = we 2 x r. (C.2) 

Since the hydrodynamic stress of a purely rotating fluid vanishes, this flow does not 
induce forces at the membrane. As the elastic energy of the membrane depends on a 
torque is exerted on the surrounding fluid. This torque is, however, of the order e 2 and 
is neglected in the following. The no-slip boundary condition at the capsule membrane 
therefore requires 

$=we 2 . (C.3) 

For v s we chose a flow which matches the elongational part of the external flow at 
infinity and vanishes at the sphere. While such a flow field produces surface forces on 
the capsule membrane, these forces only depend on the outer viscosity r\ . The resulting 
force at r = R is 

f 2 % = 5isr] (6 m ,- 2 - <W)^ (2Y 2 m + V6*™) . (C.4) 

The remaining flow field is induced by the surface forces of the capsule membrane 
and of the hydrodynamic stress jump of the extensional flow. The induced flow field 
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itself vanishes away from the capsule. In other words, the induced flow field is the same 
as in an ambient fluid. The global flow field can be found explicitly (see [89]), but is of 
no interest here. We only need to know the value of the flow field at the capsule. The 
linear relationship between the surface forces and the induced velocity is given by the 
(inverse) Oseen tensor 

The viscosity contrast only appears in the expression of the Oseen tensor, which reads 
explicitly 




^V^TT ^(2Z + 1)+77 (2Z + 1) 

Again, local area incompressibility enforces a coupling between radial and tangential 
flow. 

Finally, the membrane is advected with the resulting flow owing to the no-slip 
boundary condition. 

The membrane dynamics is thus governed by the equation 

D t u lm = v? m + O l (g m + f* m + % m ) ■ (C.7) 

Here, the material time derivative of the deformation amplitude reads 

D t u = (d t + (v • V))u. (C.8) 

In the following we only take the rotational part of the flow when calculating the 
advection, since the other contributions are smaller by a factor e. In other words, 
we use 

D t u = (d t + (* x r) ■ V)u = (d t + wfy)u. (C.9) 
This is justified because the rotational velocity is an order of magnitude larger than the 



(C.6) 



velocity components leading to a membrane deformation. Equations (B.ll, C.3, C.4 



C.5, C.7, C.9) specify the dynamics of the capsule in the quasi-spherical limit completely. 

We see that in linear external flow only the I = 2 modes can be excited. We 
therefore restrict ourselves to only taking I = 2 deformation modes. Since the radial 
and tangential modes are coupled by the local area constraint, we can eliminate the 
latter modes. By redefining the time scale r and introducing dimensionless strengths of 



the rotational A and elongational S flow analogously to [97] (see (38)) the equations of 
motions simplify to (omitting the now redundant Y and I = 2 indices) 

d T (f) =A (CIO) 

(d T - iAm)u m = (Smm> - u m u* m ,) x (Cll) 

m' 

[S^ 1 exp(-im'$)w m // sin/3 - z(8 m i )2 - 8 m >- 2 ) ■ 

These equations are formally very similar to the equations of motions for vesicles in 
(48) of |97j . Let us first search for solutions which show the same mirror symmetry with 
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respect to the xy-plane as the external flow. Thus, only the modes with m 
are non-zero (also for u m ). Since the modes are normalized to J2 m \ u m\ A 
employ the parametrisation 

no = cos(/3) 

u 2 = sin(/3)exp(-2i^)/v / 2 
U-2 = sin( / S)exp(2i*)/\/2 



-2,0,2 
1, we can 



(C.12) 
(C.13) 
(C.14) 



Here, \1/ G [— 7r/2, tt/2] is the inclination angle of the capsule with respect to the x-axis, 
(3 G [0,7r/2]. To compare with the equations of motion (21, 22) in the reduced Skotheim 
model, (cf. section [5]) it is advantageous to switch to the phase angle </> in the capsule 
frame rather than the phase angle $ in the laboratory frame. Both angles are related 
by the inclination angle, $ = + We thus arrive at the equations of motion (35 37). 
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